home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / dlae2.z / dlae2
Encoding:
Text File  |  2002-10-03  |  3.2 KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAAEEEE2222((((3333SSSS))))                                                            DDDDLLLLAAAAEEEE2222((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLAE2 - compute the eigenvalues of a 2-by-2 symmetric matrix [ A B ] [ B
  10.      C ]
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DLAE2( A, B, C, RT1, RT2 )
  14.  
  15.          DOUBLE        PRECISION A, B, C, RT1, RT2
  16.  
  17. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  18.      These routines are part of the SCSL Scientific Library and can be loaded
  19.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  20.      directs the linker to use the multi-processor version of the library.
  21.  
  22.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  23.      4 bytes (32 bits). Another version of SCSL is available in which integers
  24.      are 8 bytes (64 bits).  This version allows the user access to larger
  25.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  26.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  27.      only one of the two versions; 4-byte integer and 8-byte integer library
  28.      calls cannot be mixed.
  29.  
  30. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  31.      DLAE2 computes the eigenvalues of a 2-by-2 symmetric matrix [ A B ] [ B C
  32.      ]. On return, RT1 is the eigenvalue of larger absolute value, and RT2 is
  33.      the eigenvalue of smaller absolute value.
  34.  
  35.  
  36. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  37.      A       (input) DOUBLE PRECISION
  38.              The (1,1) element of the 2-by-2 matrix.
  39.  
  40.      B       (input) DOUBLE PRECISION
  41.              The (1,2) and (2,1) elements of the 2-by-2 matrix.
  42.  
  43.      C       (input) DOUBLE PRECISION
  44.              The (2,2) element of the 2-by-2 matrix.
  45.  
  46.      RT1     (output) DOUBLE PRECISION
  47.              The eigenvalue of larger absolute value.
  48.  
  49.      RT2     (output) DOUBLE PRECISION
  50.              The eigenvalue of smaller absolute value.
  51.  
  52. FFFFUUUURRRRTTTTHHHHEEEERRRR DDDDEEEETTTTAAAAIIIILLLLSSSS
  53.      RT1 is accurate to a few ulps barring over/underflow.
  54.  
  55.      RT2 may be inaccurate if there is massive cancellation in the determinant
  56.      A*C-B*B; higher precision or correctly rounded or correctly truncated
  57.      arithmetic would be needed to compute RT2 accurately in all cases.
  58.  
  59.      Overflow is possible only if RT1 is within a factor of 5 of overflow.
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDLLLLAAAAEEEE2222((((3333SSSS))))                                                            DDDDLLLLAAAAEEEE2222((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      Underflow is harmless if the input data is 0 or exceeds
  75.         underflow_threshold / macheps.
  76.  
  77.  
  78. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  79.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  80.  
  81.      This man page is available only online.
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.